Peierls Instability in Carbon Nanotubes 
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We present a first-principles study of Peierls distortions in irans-polyacetylene, polyacene, and 
armchair (n, n) carbon nanotubes. Our findings suggest that the ground-state geometries of arm- 
chair (n, n) carbon nanotubes, with n up to 6, exhibit a Peierls distortion as it is found for trans- 
polyactetylene. In contrast to previous studies in which no Peierls distortion is found with con- 
ventional local and semi-local density functionals, we use a hybrid functional whose exact-exchange 
admixture has been specifically optimized for the problem at hand. 



The debate whether carbon nanotubes should exhibit 
a Peierls instability arose shortly after their discovery 
Trans-polyacetylene {t-PA) is the classic example 
of Peierls instability and since (n, n) armchair nanotubes 
can be considered as a parallel arrangement of polyacety- 
lene chains along the tube axis, it is natural to suspect 
that such a deformation is also possible in these systems. 
Studies [IHil of the electron-phonon coupling using model 
Hamiltonians predict a Peierls transition with low criti- 
cal temperature. Whereas earlier investigations based on 
density-functional theory (DFT) found a symmetric ge- 
ometry, more recent investigations using denser grids or 
refined zone-folding techniques also reveal a Peierls 
instability for small nanotubes. All these studies found 
an instability at q — 2kp which results in a reduction 
of the translational symmetry. The corresponding dis- 
tortion is quite different than the deformation in f-PA 
because at least three unit cells are required to describe 
it. According to Ref. [f|, many more unit cell could be 
necessary, because they find fc^ to be incommensurate 
with the unit cell. The phonon spectrum obtained with 
DFT calculations also indicates a softening at q = 
which corresponds to a deformation that only doubles the 
unit cell, i.e. a deformation that is quite similar to the 
one found in t-PA. 

However, it has been shown 0, IH that the standard 
local density (LDA) and generalized gradient approxima- 
tions (GGAs) are unable to reproduce the experimental 
geometry of trans-polyacetylene [9l-[l2j|. underestimating 
the bond length alternation (BLA) by 70-80%. Thus, in 
order to reproduce the correct ground state of t-PA from 
a first-principles method, theories which treat exchange 
and correlation more accurately have to be employed 13- 
15j. Hybrid functionals, in which a part of the approxi- 
mate DFT exchange is replaced by exact exchange, have 
proven to give a more accurate estimat of the BLA in 
frans-polyacetyiene [Hj]. 

The previous DFT studies mentioned above used local 
and semi-local density functionals but in light of results 
for t-PA, it seems important to consider the effects of 
hybrid functionals on armchair carbon nanotubes. No- 
tably, it could change the relative strength of the insta- 



bility at q = compared to q — 2kp. Understanding 
the nature of the instability is crucial since it could well 
be related to the superconductivity found in carbon nan- 
otubes 
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In this article, we report an ab initio study of Peierls 
instabilities in irans-polyacetylene, polyacene (PAc), and 
(n, n) armchair carbon nanotubes and highlight the im- 
portance of including exact exchange in the functional to 
obtain a dimerized structure. 

Method: All calculations were done with the Kohn- 
Sham approach of DFT using gaussian basis sets, as im- 
plemented in the periodic boundary conditions (PBC) 
module of the GAUSSIAN03 package [18j . We employed 
the 6-311G and 6-31G basis sets for the polymers and 
the nanotubes respectively, and the Brillouin zones were 
sampled with over 200 k-points. Atomic positions and 
unit cell lengths of the tubes and polymers have been 
optimized. Several exchange-correlation functionals have 
been put to work in this study, namely the PBE [l9j . 
B3LYP [13, El|, HSE 0, as well as the PBE hybrid 
23]. referred thereafter as hPBE. and a generalization of 
the latter: 
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where x is the fraction of exact exchange included in 
the functional. The PBE and hPBE functionals corre- 
spond to x — and x = 0.25, respectively. We will 
refer, hereafter, to the various flavors of this functional 
by specifying only the fraction of exact exchange x. We 
would like to stress that the optimal amount of exact ex- 
change mixing depends on the system of interest 24 , 2|| ■ 



Therefore, the "best" hybrid functional for the problem 
at hand contains one unknown parameter. Furthermore 
we should mention that the implementation of the hy- 
brid methods in the GAUSSIAN does not strictly adhere 
to the Kohn-Sham approach since the non-local Hartree- 
Fock exchange operator is employed instead of the local 
exchange potential. Experience indicates that this modi- 
fication has essentially no effect on the predicted ground- 
state properties [13]. The HSE functional used in this 
study includes 25% of i^ xa ct an( j ^ e Coulomb interac- 
tion is screened at a distance of 3.5 A (see the reference 
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FIG. 1. Most stable configurations obtained for PAc and 
(n, n) CNTs with x > 0. Here the double, aromatic and sin- 
gle bonds notation is used to emphasize that the bond lengths 
are ordered as follows A < B < C. 

by Heyd et al. [22J for details of the implementation) . 

Distortion vs exact exchange: In this section we re- 
port the bond length alternation along the system axis, 
defined as the difference between the long B and short 
A bonds, and the associated band gap obtained using 
different exchange-correlation functionals. The systems 
considered are frons-polyacetylene, polyacene, and car- 
bon nanotubes. The distortion found to be the most sta- 
ble in PAc and (n, n) tubes is depicted in Fig. [TJ These 
systems are very similar to that of a parallel arrange- 
ment of distorted i-PA chains along the tube (or poly- 
mer) axis. The BLA of neighboring chains exhibits a 
phase difference corresponding to half a cell. Hence, the 
distorted geometry is found to be a Kekule-type struc- 
ture, the bond lengths being all different and are ordered 
as follows: A < B < C (see Table HH). Note that it is 
also different from what was reported in [la], but it is very 
similar to that of [l[ and the one found in [27| , which also 
found a Kekule-type structure. 

TableU shows the BLA (denoted by Aa) and the band 
gap, E g , obtained with the B3LYP (20% of £f xact ), 
hPBE (25% of £| xact ) and HSE (25% of screened 
^exact) f unc tionals for t-PA, PAc and the (3,3) nan- 
otube. As expected, the band gap is proportional to the 
amount of exact exchange included in the functional and 
so is the BLA. This stems from the fact that Hartree- 
Fock theory tends to overestimate band gaps and over 
localize electrons. Hence, an appropriate choice of the 
parameter x can recover the correct BLA. We see from 
Table U that although the HSE functional yields a dimer- 
ized structure for t-PA and PAc, the magnitude of the 
BLA is less than that obtained with functionals which in- 
clude the full extent of the exchange interaction suggest- 
ing that the medium to long-range exchange interaction 
increases the dimerization in these systems. In particu- 
lar, for the (3,3) nanotube, HSE predicts a symmetric 
geometry with a very small gap. The importance of this 
long-range interaction could explain the failure of purely 
local- and semi-local functionals, such as LDA and PBE. 

To systematically study the exact exchange depen- 
dence of the dimerization amplitude and band gap, we 
used the functional of Eq. [1] and varied x between 





B3LYP 


hPBE 


HSE 




Aa Eg 


Aa 


E g 


Aa E g 




(10" 3 A) (eV) 


(io- 3 A) 


(eV) 


(10" 3 A) (eV) 


f-PA 


56.4 1.23 


59.0 


1.49 


49.6 0.84 


PAc 


24.8 0.44 


26.7 


0.55 


17.0 0.17 


(3,3) 


4.13 0.25 


4.97 


0.36 


0.0 0.05 



TABLE I. Calculated distortions and band gaps for t-PA, 
PAc, and a (3,3) tube with different exchange-correlation 
functionals. 

and 0.5. Fig. [5] shows the dimerization amplitude (top 
panel) and the band gap (bottom panel) of t-PA, PAc, 
and (n,n) CNTs, for n ranging from 3 to 5, as a func- 
tion of x. The geometries were optimized for each value 
of x. The reported band gaps are differences in Kohn- 
Sham eigenvalues and should be treated only as indica- 
tions of a trend since the electron-hole interaction is not 
included and it has been shown to be important in nan- 
otubes [2^, [29[. The BLA in the tubes scales linearly 
with the amount of exact exchange, while the band gap 
follows a power law x m with m > 1. This figure also 
shows that the density only PBE (x = 0) functional is 
unable to reproduce the experimental BLA of 0.08 A nor 
the band gap of 1.5 eV [30j, l3l( of frans-polyacetylene 
and predicts all (n, n) CNTs to be conductors. The in- 
crease in the band gap as a function x is caused by the 
addition of exact exchange included in the functional and 
by the geometry of the tube that is more dimerized as x 
increased. 

Optimized geometry of the tubes: To study the Peierls 
mechanism in more detail, we use the amount of ex- 
act exchange that best describes the frcms-polyacetylene 
ground state geometry since the tubes and PAc can be 
viewed as an arrangement of f-PA chains. For PAc the 
best agreement with the experimental geometry was ob- 
tained for x = 0.3. We prefer this hybrid functional over 
the B3LYP (or a generalization of it) because it recov- 
ers the homogenous electron gas limit. In Table ITT1 the 
optimized bond lengths, differences in total energy with 
respect to the symmetrical case and the band gaps are 
reported for t-PA, PAc and (n, n) obtained with x = 0.3. 
From these results it is clear that the distortion becomes 
less significant as the tube diameter increases. Also, all 
bond lengths seem to converge to the same value of 1.421 
A which corresponds to the hPBE(a; = 0.3) value for the 
C-C bond length in graphite. 

The stability of these deformations was evaluated by 
calculating the total energy gain achieved by deform- 
ing the tube. The undistorted geometry was chosen 
such that it has equal A and B bond lengths, a com- 
mon feature of aromatic systems. Note that within LDA 
and PBE, armchair carbon nanotubes also exhibit equal 
bond lengths. The structure was permitted to relax, 
under the above mention condition (A=B). This undis- 
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FIG. 2. (color online) Dimerization amplitude (top panel) 
and band gap (bottom panel) as a function of the percentage 
of exact exchange (a;) included in the functional. The shaded 
area indicates the range of x yielding the best agreement be- 
tween calculated and experimental £-PA bond lengths. The 
insets zoom in on the small- a; region. 





Bond lengths (A) 




n 


A 


B 


C 


A^tot 


(meV, K) 


E£< B (eV) 


f-PA 


1.361 


1.426 




27.11 


314.56 


1.66 


PAc 


1.383 


1.416 


1.460 


15.46 


179.44 


0.749 


3 


1.426 


1.432 


1.440 


3.41 


39.57 


0.49 


4 


1.423 


1.426 


1.431 


2.10 


24.44 


0.36 


5 


1.422 


1.424 


1.427 


1.14 


13.22 


0.27 


6 


1.421 


1.423 


1.424 


0.69 


7.98 


0.23 



TABLE II. Optimized geometries of t-PA, PAc and (n, n) 
CNTs with x = 0.3. 



torted structure reaches its lowest total energy state when 
= B(0) = mean(A(P er \ J B(P er )). The bond-length C 
was allowed to relax but it was found to remain the same 
as in the fully-optimized structure. We calculated the to- 
tal energy difference, AE^ ^, between the distorted and 
undistorted structures for all tubes. We report the re- 
sults of this analysis in Table [TTJ The stability of the 
deformations decreases rapidly with increasing tube di- 
ameter. Using these energy differences as an estimation 
of the temperature Tpo at which the PD would occur, 
we find that it corresponds to a few kelvins for the nan- 
otubes. Our estimates are between what was predicted 



by Bohnen et al. Q and Connetable et al. [B[ , i.e. a near 
room temperature transition for the (3,3) nanotube, and 
what was suggested by Mintmire et al. [lj and Piscanec 
et al. 0, i.e. T PD < 1 K. 

To our knowledge, the Peierls distortion in carbon nan- 
otubes has never been directly and conclusively observed 
in experimental data, owing much from the fact that most 
readily available nanotubes are larger then the one stud- 
ied here. As we have shown, the stability of such dis- 
tortions decreases rapidly with increasing tube diameter. 
As such, it is comprehensible that small perturbations 
are enough to wash out this effect. Nonetheless, our re- 
sults seems to concur with the resistance measurements 
of Shea et al. [32| , which found a sharp increase in resis- 
tance at T = 0.8K for a bundle of tubes. The diameter 
of the most abundant nanotube in the sample was 1.4 
nm, corresponding to a (10,10) tube. Extrapolating our 
results to such a tube, using a simple exponential fit [33j . 
we find T PD = 0.9 K. 

Connection with the response function: Within linear 
response theory, the Peierls instability can be understood 
as a pole in the response function \(r, r') or equivalently 
as a zero eigenvalue in the X i r -> r') spectrum. In a DFT 
approach, the real response of the system is related to 
the response function of the non-interacting Kohn-Sham 
system denoted by X s { r i r ') by the following relation: 



_J_ = 1 _ + f U r i)\ , 

X(r,r') X s (r,r') \\r-r'\ ' / ' 



(2) 



where / lc (r,r') is the exchange-correlation kernel. For a 
one-dimensional system, the Kohn-Sham non-interacting 
electron system exhibit a Peierl instability at a wave vec- 
tor q = 2kp. The Coulomb part of the previous equation 
has a negative contribution. Since the response function 
X s {r, r') is negative definite, this Coulomb part will make 
the spectrum of x( r i r ') more negative. Hence it tends to 
decrease any instabilities in Kohn-Sham non-interacting 
electron system. On the other hand, the exchange- 
correlation kernel, since it is always negative, should at- 
tenuate/cancel the Coulomb part and preserve the Peierls 
instability. Therefore, the exchange-correlation kernel 
should contain a long-range contribution of the form — 
where a is a material dependent parameter. We note that 
a kernel of the form x(q) = has been suggested for 
TDDFT calculations in order to obtain correct optical 
spectra [34j. Here we find that such response functions 
are important for the ground-state geometry in the sys- 
tem studied in this paper. 

Based on this discussion, it is not surprising that ex- 
act exchange helps to reestablish the Peierls instability. 
Essentially, the direct Coulomb interaction alone would 
diminish the Peierls instability since it would tend to re- 
frain electrons to localize in a same region. The local 
and near local functionals (LDA, PBE and HSE) tend 
to a constant has q — > 0. Hence, they fail to produce 
an appreciable Peierls distortion as they cannot cancel 
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the Coulomb contribution. The inclusion of some exact- 
exchange restores a correct long- wave dependence of the 
kernel. The amount of exact-exchange must be chosen as 
to give the correct a parameter, which unfortunately will 
depend on the system studied. In this paper, we suggest 
that 30% of exact- exchange is nominal for the study of 
ID conjugated carbon systems such as t-PA, PAc and 
armchair nanotubcs. 

Conclusion: Our DFT calculations show that hybrid 
functionals are better suited for the description of the 
properties of trans-polyacetylene than local- or semi-local 
functionals (LDA or GGAs), and hence could be more ap- 
propriate to describe the geometry of carbon nanotubes. 
Using a hybrid functional derived from PBE, we showed 
that the distortion amplitude and the band gap increases 
with the amount of exact exchange included in the func- 
tional. Using the value of x that best describes the t- 
PA geometry, we found that small radius armchair nan- 
otubes should exhibit small lattice distortions that de- 
crease with increasing tube diameter. Our particular 
choice of functional appears appropriate since it recovers 
the graphite bond lengths for larger tubes, in addition to 
describing the bond length alternation in t-PA correctly. 
We estimated the Peierls transition temperature to be of 
the order of few degrees K. 
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